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Abstract 

A search for scalar and vector leptoquarks coupling to first generation fermions is per- 
formed using the e + p and e~p scattering data collected by the HI experiment between 
1994 and 2000. The data correspond to a total integrated luminosity of 117pb _1 . No evi- 
dence for the direct or indirect production of such particles is found in data samples with a 
large transverse momentum final state electron or with large missing transverse momentum. 
Constraints on leptoquark models are established. For leptoquark couplings of electromag- 
netic strength, leptoquarks with masses up to 275 — 325 GeV are ruled out. These limits 
improve and supercede earlier HI limits based on subsamples of the data used here. 
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Introduction 



The ep collider HERA offers the unique possibility to search for the resonant production of 
new particles which couple directly to a lepton and a parton. Leptoquarks (LQs), colour triplet 
bosons which appear naturally in various unifying theories beyond the Standard Model (SM), 
are such an example. At HERA, leptoquarks could be singly produced by the fusion of the 
initial state lepton of energy 27.6 GeV with a quark from the incoming proton of energy up to 
920 GeV. 

The phenomenology of LQs at HERA is discussed in detail in [1]. The effective Lagrangian 
considered conserves lepton and baryon number, obeys the symmetries of the SM gauge groups 
U(l)y, SU(2)l and SU(3)c and includes both scalar and vector LQs. A dimensionless param- 
eter A defines the coupling at the lepton-quark-LQ vertex. At HERA, LQs can be resonantly 
produced in the s-channel or exchanged in the u-channel between the incoming lepton and a 
quark coming from the proton. In the s-channel, a LQ is produced at a mass M = ^Jxs ep where 
x is the momentum fraction of the proton carried by the interacting quark. 

This letter presents a search for LQs coupling to first generation fermions using e + p data col- 
lected at a centre-of-mass energy of yfs^, « 300 GeV, e~p data collected at v /s^ ~ 320 GeV 
and e + p data collected at ^fs^, « 320 GeV. These data sets correspond to integrated luminosi- 
ties of 37 pb" 1 , 15 pb _1 and 65 pb" 1 respectively. They represent the full statistics accumulated 
by the HI experiment between 1994 and 2000. The results of this search thus supercede those 
based on the e+p 1994-1997 [1] and e~p 1998-1999 [2] data. 

The e + p and e~p data are largely complementary when searching for leptoquark resonances. 
Due to the more favourable parton densities of quarks with respect to anti-quarks at high x, the 
e + p (e~p) data sets are mostly sensitive to LQs with fermion number 1 F — (F — 2). The 
search reported here considers the leptoquark decays into an electron and a quark and into a 
neutrino and a quark. These LQ decays lead to final states similar to those of deep-inelastic 
scattering (DIS) neutral current (NC) and charged current (CC) interactions at very high Q 2 , the 
negative four-momentum transfer squared. If the final state consists of an electron and a quark, 
the LQ mass is reconstructed from the measured kinematics of the scattered electron. If the 
final state consists of a neutrino and a quark, the LQ mass is reconstructed from the hadronic 
final state [1]. 



Neutral and Charged Current Data 

The HI detector components most relevant to this analysis are the liquid argon calorimeter, 
which measures the positions and energies of charged and neutral particles over the range 
4° < 9 < 154° of polar angle 2 , and the inner tracking detectors which measure the angles 
and momenta of charged particles over the range 7° < 9 < 165°. A full description of the 
detector can be found in [3]. 

'The fermion number F is given by F = 3B + L with B and L being the baryon and lepton numbers respec- 
tively. 

2 The polar angle 9 is defined with respect to the incident proton momentum vector (the positive z axis). 
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This search is based on inclusive NC and CC DIS data in the kinematic domain Q 2 > 
2500 GeV 2 and 0.1 < y < 0.9, where the inelasticity variable y is defined as y = Q 2 /M 2 . The 
cuts on y remove regions of poor reconstruction, poor resolution, large QED radiative effects 
and background from photoproduction processes. The selection of NC-like events follows that 
presented in [1]. It requires an identified electron with transverse momentum above 15 GeV. 
The selection of CC-like events follows closely that presented in [1, 4]. A missing transverse 
momentum exceeding 25 GeV is required. 

The inelasticity variable y is related to the polar angle 9* of the decay lepton in the centre- 
of-mass frame of the hard subprocess (eg — > Iq) by y = |(1 + cosO*). Since the angular 
distribution of the electron coming from the decay of a scalar resonance is markedly different 
from that of the scattered lepton in NC DIS [1], a mass dependent cut on y was applied previ- 
ously [1, 2] in order to optimise the signal sensitivity. However, the optimisation power is rather 
limited for a vector resonance as the angular distribution is only slightly different from that of 
the DIS background. For this reason, no such mass dependent y cut is applied in this analy- 
sis. Instead, all selected events are analysed in bins of varying size adapted to the experimental 
resolution in the M — y plane, with a procedure designed to fully exploit the sensitivity to the 
signal as explained in the next section. 

The mass spectra measured for NC- and CC-like events in the three data sets are compared in 
Figs, la-f with the SM predictions, obtained using a Monte-Carlo (MC) calculation [5] and the 
CTEQ5D parametrisation [6] for the parton densities. In all cases the data are well described 
by the SM prediction. Since no evidence for LQ production is observed in either the NC or 
CC data samples, the data are used to set constraints on LQs which couple to first generation 
fermions. 



For the limit analysis, the data are studied in bins in the M — y plane. The binning used is 
different for the different data sets and is shown in Fig. 2. For those LQ types with only a NC- 
like decay mode, the total number of bins amounts to about 200 covering all three NC samples. 
The total number of bins doubles when including the three CC data samples for LQs having 
both NC-like and CC-like decay modes. The number of SM background events bi in each bin 
i is obtained from the SM MC calculations. Each MC event k, reconstructed in bin i, has an 
event weight e^, such that the MC is normalised to the luminosity of the data. The sum over all 
SM MC events within bin i thus gives 



To estimate the LQ signal, an event re-weighting technique is applied to the same SM MC 
events. No use is made of a dedicated signal MC generator. The number of events expected in 
each bin i in the presence of a LQ signal is denoted as s, + b, L . It may be written as 



Statistical Method 




(i) 




(2) 
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where a -£Q> INT > SM are differential cross section terms [7] corresponding to the LQ, interference 
and SM contributions, respectively. These differential cross section terms, calculated in leading 
order using the parton density functions CTEQ5D [6], are based on the true kinematic quan- 
tities of event k, whereas the resulting Sj + bi events are counted in the M — y bins of the 
reconstructed variables with appropriate simulation of the detector response. The differential 

LO INT 

cross section terms a k depend on the LQ mass and coupling A. For mass values well be- 
low the kinematic limit yfs^,, the s-channel contribution dominates in the LQ cross section and 
the signal contribution s, to bin i is always positive. However, at higher masses, the interfer- 
ence contributions become more important and Si may be negative in the case of destructive 
interferences, although s, + bi always stays positive. 

The limits are determined from a statistical analysis which uses the method of fractional 
event counting [8]. For a given leptoquark mass and coupling a weight Wi is ascribed to each 
bin, which is given by the asymmetry between the expected number of events in the presence 
or absence of a LQ signal: 

w _ (sj + bi) - h _ Si 

(si + k) + bi Si + 2bi ' 

As an example, Fig. 2 illustrates the bin dependent weights for an e + d type vector LQ having 
a mass of 200 GeV, a coupling of 0.023 and both the NC-like and CC-like decay channels. As 
expected, the weights are small for the e~p data. For the e + p data sets, the weights have little 
y-dependence. 

Using these weights a fractional event count, also called test statistics, is defined as 

X(data) = Y^iV^data) , (4) 

i 

where JVj(data) is the number of data events observed in bin i. 

In a frequentist approach, a large number of "experiments" (2 x 10 000) are generated. Each 
experiment consists of Poisson distributed random numbers Ni(b) (Ni(s + b)), based on the 
expected number of events bi (s» + bi) in the absence (presence) of a LQ signal. For each 
background experiment b and for each signal-plus-background experiment s + b a fractional 
event count is defined in analogy to Eqn.(4): 

X(b) = ^WiNiib) (5) 

i 

X( S + b) = J2 W i N i( S + b )- ( 6 ) 

i 

Frequentist probabilities CL s+b (CL&) are defined as the fraction of experiments where the quan- 
tity X(s + b) (X(b)) is smaller than X(data). If the data agreed perfectly with the expectation 
from the background-only hypothesis, a value of CL & = 0.5 would be obtained. A higher value 
indicates that the observation is more signal-like; a lower value indicates fluctuations opposite 
to those expected for a signal. If CL S+ ;, is small, it may be used to exclude the signal-plus- 
background hypothesis with confidence level (1 — CL S+ ;,). However, in this analysis we use the 
confidence level CL defined as 

CL = 1 - CL s+b /CL b (7) 
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to set limits in a conservative manner. This ratio, which was also used in LEP searches [8, 9], has 
the desirable feature that as the LQ coupling tends to zero, and necessarily CL s+fe — > CL fe , CL 
drops to 0, i.e. one cannot rule out any LQ which has a vanishing coupling. On the other hand, 
for a non-zero coupling, an exclusion limit at 95% CL is always reachable for a sufficiently 
large coupling. 

Systematic uncertainties enter as offsets 6^ and 5^ b to the predicted number of events 6j 
and Si + bi, where j runs over all independent sources of systematic errors. For the limits 
presented below, both the bin weights Wi and the b and s + b MC experiments are altered by the 
known systematic uncertainties, assuming that the latter have Gaussian probability densities. 

The experimental systematic error is dominated by the electromagnetic energy scale (be- 
tween 0.7% and 3%) for the NC analysis, and by the hadronic energy scale (2%) for the CC 
analysis. The limited knowledge of proton structure causes an uncertainty on the signal cross 
section. This uncertainty is estimated to be 5% for F = 2 (F = 0) LQs coupling to e~u (e + u) 
and varies between 7% at low LQ masses up to 30% around 290 GeV for F = 2 (F = 0) LQs 
coupling to e~d (e + d). Similarly, an uncertainty on the DIS cross sections is connected with 
the parton densities. The correlation between the systematic uncertainties on the signal and that 
on the background and between different analysis bins is taken into account. The correlations 
induced by the uncertainties of the parton densities are evaluated using [10]. 



Limit Results 

In the following limits will first be derived within the phenomenological model proposed by 
Buchmuller, Ruckl and Wyler (BRW) [7] and then within generic models where the branching 
ratios (3 e (J3 V ) for the LQ decays into eq (uq) are not fixed. 

The BRW model describes 7 LQs with F = and 7 LQs with F — 2. We use here the 
nomenclature of [1 1] to label the various scalar SV,l OSj r) or vector Vj L (V ItR ) LQ types of 
weak isospin I, which couple to a left-handed (right-handed) electron. The tilde is used to 
distinguish LQs which differ only by their hypercharge. In the BRW model the branching ratios 
P e (J3 V ) are fixed and equal to 1 or 0.5 (0 or 0.5) depending on the LQ quantum numbers. Table 1 
lists the 14 LQ types according to the BRW model. 

For LQs with F — 0, the upper limits on the coupling obtained at 95% CL are shown as a 
function of the LQ mass in Figs. 3a and b, for scalar and vector LQs respectively. For masses 
above ~ 270 GeV, these bounds improve by a factor of up to ~ 3 the limits obtained in [1] 
from the analysis of e + p data at y/s^ = 300 GeV. Constraints corresponding to F = 2 LQs are 
shown in Figs. 3c and d which extend those in [2] beyond the kinematic limit. For a coupling 
of electromagnetic strength a cm (A = \JAiia crn = 0.3) this analysis rules out LQ masses below 
275 to 325 GeV, depending on the LQ type. 

Fig. 4 summarises the constraints on the Su2,l and on the So,l obtained by HI, by the OPAL 
and L3 experiments at LEP [12], and by the DO experiment at the Tevatron [13]. The limits 
shown from LEP are from indirect constraints. The limits from the Tevatron are independent 
of the coupling A as they were derived from the dominant pair production processes. The HI 
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Table 1 : Leptoquark isospin families in the Buchmiiller-Riickl-Wyler model. Charge conjugate pro- 
cesses are not shown. For each leptoquark, the subscript denotes its weak isospin. For simplicity, the 
leptoquarks are conventionally indexed with the chirality of the incoming electron which could allow 
their production in e~p. The variable (3 e denotes the branching ratio of the LQ into e + q. 



limits are comparable with those from the ZEUS experiment obtained in a similar analysis [14]. 
The limits at high mass values are also compared with those obtained in a contact interaction 
analysis [15], which was based on the measured single differential cross sections da/dQ 2 from 
the NC process only [16]. 

Beyond the BRW ansatz, generic LQ models can also be considered. An example is pro- 
vided by supersymmetric models 3 where the R-parity is violated by a Ay X (A' llfc ) coupling, with 

the u 3 L (d R ) squark having the same interactions with a lepton-quark pair as the Si/ 2 ,l (So,l)- 
In generic LQ models other LQ decay modes are allowed such that the branching ratios (3 e and 
(5 V are free parameters. Mass dependent constraints on the LQ branching ratios can then be set 
for a given value of A. For a vector LQ coupling to e + d (possessing the quantum numbers of 
the Vq : l) an d for A = 0.06, a domain of the (3 e -M ((3 V -M) plane can be excluded by the NC 
(CC) analysis as shown in Fig. 5a. If the LQ decays into eq or vq only 4 , the combination of 
both channels rules out the part of the plane on the left of the second full curve from the left 
for A = 0.06. The resulting combined bound is largely independent of the individual values 

3 More general limits on squark production taking direct and indirect i?-parity violating decay modes into 
account have been set in [17]. 

4 It should be noted that j3 e + /3„ = 1 does not imply (3 e = /3 V even when invariance under SU{2)l transfor- 
mations is required. For example, when LQs belonging to a given isospin multiplet are not mass eigenstates, their 
mixing usually leads to different branching ratios in both channels for the physical LQ states. 
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of (3 e and (3 V . Combined bounds are also shown for A = 0.03 and A = 0.3. For a coupling 
A = 0.3 and high (3 V the limit extends to high mass values above the kinematic limit of resonant 
LQ production. For this part of the parameter space, the coupling \ u = X^/f3 u /f3 e is large 5 but 
still satisfies A^/47r < 1. A smooth transition is observed between limits driven by resonant 
production and limits driven by contact interactions. Fig. 5b shows similar exclusion limits as 
for Fig. 5a, for a scalar LQ possessing the quantum numbers of the S 0j l (which couples to e~u). 
The domain excluded by the DO experiment at the Tevatron [13] is also shown. For A greater 
than ~ 0.06, the HI limits on scalar LQs extend considerably beyond the region excluded by 
the DO experiment. 

To summarise, a search for leptoquarks with fermion numbers F = 2 and F = has been 
performed using all e + p and e~p data collected by HI between 1994 and 2000. No signal has 
been observed and constraints on leptoquarks have been set, which extend beyond the domains 
excluded by other experiments at LEP and the Tevatron. For a coupling of electromagnetic 
strength, leptoquark masses below 275 — 325 GeV, depending on the leptoquark type, can be 
ruled out. 
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Figure 1 : Mass spectra for the (a-c) neutral current (NC) and (d-f) charged current (CC) deep 
inelastic scattering selected events, together with the corresponding Standard Model (SM) ex- 
pectations. The shaded bands indicate the ±ler uncertainty on the SM expectations. 
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Figure 2: Binning used in the M — y plane for the different data sets and weights calculated 
in these bins for a 200 GeV vector leptoquark with a coupling of 0.023 to e + d and vu. The 
left plots correspond to the neutral current (NC)-like decay channel whereas the right plots 
correspond to the charged current (CC)-like decay channel. The top, middle and bottom plots 
correspond to the 300 GeV e + p, 320 GeV e~p and 320 GeV e + p data sets, respectively. 
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Figure 3: Exclusion limits for the 14 leptoquarks (LQs) described by the Buchmuller, Ruckl 
and Wyler (BRW) model. The limits are expressed at 95% CL on the coupling A as a function 
of the leptoquark mass for the (a) scalar LQs with F — 0, (b) vector LQs with F — 0, (c) scalar 
LQs with F = 2 and (d) vector LQs with F — 2. Domains above the curves are excluded. 
For each LQ type the pairs of Standard Model fermions coupling to it are indicated in brackets 
(charge conjugate states are not shown). 
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Figure 4: Exclusion limits at 95% CL on the coupling A as a function of the leptoquark (LQ) 
mass for S\/2,l (top) and Sq,l (bottom) in the framework of the BRW model. The direct DO 
limits are independent of the coupling. For Sx/ 2 ,l the indirect limit from OPAL is shown, 
whereas for S 0j l the better indirect limit from L3 is shown. Constraints on LQs with masses 
above 350 GeV obtained from the HI contact interaction (HI CI) analysis [15] are also shown, 
in the rightmost part of the figures. 
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Figure 5: (a) Domains ruled out by the combination of the NC and CC analyses, for a vector 
LQ which couples to e + d (with the quantum numbers of the Vo l) and decaying only into eq 
and vq for three values of the coupling A. (b) Same as for (a) but for a scalar LQ coupling 
to e~u (with the quantum numbers of the So,l)- The regions on the left of the full curves are 
excluded at 95% CL. For A = 0.06, the part of the j3 e - M LQ (j3 v - M LQ ) plane on the left 
of the dashed (dotted) curve is excluded by the NC (CC) analysis alone. The branching ratios 
(3 e and /?„(= 1 — j3 e ) are shown on the left and right axes respectively. The excluded domains 
cover p e values larger than 7.2 x 10" 5 , 2.9 x 10~ 4 and 7.1 x 10" 3 for A = 0.03, 0.06 and 0.3 
respectively. In (b) the hatched region represents the domain excluded by the DO experiment. 
The DO limits do not depend on the value of the coupling. 
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